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Abstract. - Bohmian mechanics is a causal interpretation of quantum mechanics in which 
particles describe trajectories guided by the wave function. The dynamics in the vicinity of 
nodes of the wave function, usually called vortices, is regular if they are at rest. However, 
vortices generically move during time evolution of the system. We show that this movement 
is the origin of chaotic behavior of quantum trajectories. As an example, our general result is 
illustrated numerically in the two-dimensional isotropic harmonic oscillator. 



De Broglie-Bohm's (BB) approach to quantum mechanics has experienced an increased 
popularity in recent years. This is due to the fact that it combines the accuracy of the 
standard quantum description with the intuitive insight derived from the causal trajectory 
formalism, thus providing a powerful theoretical tool to understand the physical mechanisms 
underlying microscopic phenomena [1,2]. Although the behavior of quantum trajectories 
is very different from classical solutions it can be used to gain intuition in many physical 
phenomena. Numerous examples can be found in different areas of research. In particular, we 
can mention studies of barrier tunneling in smooth potentials [3] , the quantum back-reaction 
problem [4] and ballistic transport of electrons in nanowires [5] . 

According to the BB theory of quantum motion, a particle moves in a deterministic orbit 
under the influence of the external potential and a quantum potential generated by the wave 
function. This quantum potential can be very intricate because it encodes information on 
wave interferences. Based on it, Bohm already predicted complex behavior of the quantum 
trajectories in his seminal work [6]. This was recently confirmed in several studies when pres- 
ence of chaos in various systems has been shown numerically [7-9]. However, the mechanisms 
that cause such a complex behavior is still lacking. In this letter we show that movement of 
the zeros of the wave function, commonly known as vortices, implies chaos in the dynamics of 
quantum trajectories. Such a movement perturbs the velocity field producing transverse ho- 
moclinic orbits that generate the well known Smale horseshoes which is the origin of complex 
behavior. Our assertion is based on an analytical proof in a simplified model which resem- 
bled the velocity field near the vortices. In addition, we present a numerical study in a 2-D 
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isotropic harmonic oscillator that displays a route to chaos dominated by this mechanism. It 
is important to mention that there is no agreement in previous works about the influence of 
vortices on the chaotic motion of quantum trajectories [8-11]. 

The fundamental equations in the BB theory are derived from the introduction of the wave 
function in polar form, ip(r,t) — R(r,t) e lS ^ r,t ^ (throughout the paper h is set equal to 1), into 
the time-dependent Schrodinger equation, thus obtaining two real equations: 

»* W*^-0, (1) 



dt 



rn 



which are the continuity and quantum Hamilton-Jacobi equations, respectively. The last 
term in the left-hand-side of Eq. (J2J) is the so-called quantum potential, a non-local function 
determined by the quantum state, which, together with V, determines the total force acting 
on the system. Similarly to what happens in the usual classical Hamilton-Jacobi theory, 
quantum trajectories of a particle of mass m can then be defined by means of the following 
velocity field equation: 

V = r= m V5= 2^ W ' () 

Vortices appear naturally in the BB framework. They result from wave function inter- 
ferences so they have no classical explanation. In systems without magnetic field, the bulk 
vorticity V x v in the probability fluid is determined by the points where the phase S is singu- 
lar. This may occur only at points where the wave function vanishes. This condition is fulfilled 
by isolated points in a 2-D system and lines in a 3-D system. Due to the single- valuedness of 
the wave function, the circulation T along any closed contour £ encircling a vortex must be 
quantized, that is, 

T = / rdr = , (4) 



with n an integer [12,13]. So, the velocity v must diverge as one approaches a vortex. In fact, 
the time dependent velocity field in the vicinity of a vortex located at time t in r v (t) is given 
by 

r v (t) x w x w* 

(5) 



2l> ' \r-r v (t) 1 



.w 



where w = Vip(r v (t)) [10, 14]. We consider here 2-D systems but we belive that our results 
are valid for systems of higher dimensions. 

Before presenting our analytical results we show a numerical simulation of quantum trajec- 
tories in a system consisting of a particle of unit mass in a 2-D isotropic harmonic oscillator. 
We have set the angular frequency u> = 1, so the Hamiltonian of the system results 

1 , 8 2 d 2 1 



H =-¥^ + w )+ ^ 2+y2) - (6) 

The eigenenergies are E nxrly = n x + n y + 1 and the eigenfunctions (j> n:c n y (x, y) = exp(— |(x 2 + 
V 2 )) H nA x ) H n y {y) / \Jk 2 n * +n y n x \n y \ with n x = 0, 1, n y = 0, 1, ... where H n is the n-th 
degree Hermite polynomial. 
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Fig. 1 - Poincare surface of section for the quantum trajectories generated by the wave function of 
Eq. with b = c and a/b = (a) , 0.0553 (b), 0.1138 (c) and 0.17651 (d) with fixed value of 
71 = 3.876968 and 72 = 2.684916. The trajectories of the vortex of the corresponding wave functions 
are shown in Fig. |3 

Fig. 2 - Part of the Poincare surface of section for the quantum trajectories generated by the wave 
function of Eq. J7J with small values of a/b. The parameter b = c and a/b = 0.01082 (a) , 0.02175 
(b), 0.0328 (c) and 0.0440 (d) with fixed value of 71 = 3.876968 and 72 = 2.684916. 



We have chosen the following general combination of the first three eigenstates of the 
Hamiltonian of Eq. © as initial state 

ipo = a^oo + 6 exp(— 171)^10 + cexp(— 172)^01, (7) 

with a, 6, c, 71 and 72 real numbers and a 2 + b 2 + c 2 = 1 (the normalization condition). A 
remarkable point is that this state generates a periodic time dependent velocity field with only 
one vortex. Moreover, the trajectory of the vortex can be obtained analytically resulting 

... , a sin(7 2 — t) a sin(7i — t) „ 

r v (t = (x v (t ),y v (t = — ^ '-, . , U1 '-). (8) 

V26 sm(7i - 72 ) V2c sin(7i - 72 ) 

This fact allows us to see the influence of the movement of a vortex in the dynamics of the 
quantum trajectories, without taking into account the possibility of instantaneous creation or 
annihilation of a vortex pair with opposite circulation [14,15]. This important phenomenon 
will be studied elsewhere [16]. 

The non-autonomous velocity field generated by the wave function of Eq. (JJJ is periodic 
so the best surface of section is given by fixing t — 2nn with n — 0,1,.... (also called a 
stroboscopic view). Fig^shows surfaces of section with b = c and a/b = 0, 0.0553, 0.1138 and 
0.17651 with fixed value of 71 = 3.876968 and 72 = 2.684916. The trajectories of the vortex 
for the cases studied in Fig. are plotted in Fig. [3] A clear transition to chaos appears as 
the parameter a/b is increased. If the position of vortex is fixed the trajectories are regular 
and no chaos is present [see Fig.^a)] . However, irregular dynamics is observed for small a/b 
[Fig. njb)]. The transition to irregular dynamics is shown in Fig. El The movement of the 
vortex produces a saddle point near (0.6,0.75) and their stable and unstable manifolds have 
a topological transverse intersection generating the well known homoclinic tangle [17]. 
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Fig. 3 - Path described by the vortex [Eq. JHJ] °f the velocity field generated by wave functions of 
Eq. with b = c and a/b = (filled circle) , 0.0553 (dotted line), 0.1138 (dashed line) and 0.17651 
(solid line) with fixed value of 71 = 3.876968 and 72 = 2.684916. 

Now we will show analytical results that explain the numerical experiments presented 
before. Our starting point is the following model: a particle of unit mass on the plane in the 
velocity field of Eq. JSJ with the constraints that the trajectory of the vortex is a time periodic 
curve and w x = iw y [18]. Thus, the non-autonomous vector field is equal to 



-(y-yv(t)) 



(x-x v (t)f + {y-y v {t)f 
(x - x v (t)) 



y {x~x v {t)f + {y-y y {t)f 

Taking x — x — x v (t)y = y — y v (t) and writing in polar coordinates (x = rcos(#), y 
■sm(6)), the velocity field of Eq. results 



(9) 



v r = r[sin(9)y v (t) + x v (t) cos(9)] 

v e = - + cos(9)y v (t) - x v (t) sin(6>). 
r 

This non-autonomous velocity field can be seen as a perturbation of the autonomous veloc- 
ity field v ee (0, i), with G(r, 0, t) = (r[sm{0)y v (t) + x v {t) cos(9)},cos(0)y v (t) - x v (t) sm(9)) 
the time-periodic perturbation. Note that the field is induced by the time dependent Hamil- 
tonian 



H{r,0,t) = -log(r)+r[cos(e)y v (t)-sm(6)x v (t)]. 



(10) 



We consider periodic curves r v (t) such that 



cos(9)y v (t) — x v (t) sin(9)dtds ^ 0. 



(11) 



Under these hypothesis the following property can be proved: 

There exists a saddle periodic orbit of the flow associated to the vector field of Eq. 0). 
exhibiting a homoclinic transversal intersection. 

This result, which implies that quantum trajectories show topological chaos, is the main 
finding of our work . We illustrate here the geometrical arguments of the proof and we leave 
the full details for a future publication [19]. We notice that a periodic orbit of saddle type 
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Fig. 4 - (a) Schematic plot of the invariants of the non-perturbed map R° generated by the autonomous 
velocity field vo (solid lines). With dashed line it is plotted a segment with 8 = 0. (b) Non-perturbed 
mapped R° of the segment with 9 = (dashed line), (c) It is showed with dashed line the image of 
the perturbed map 7? of a circle CV with radius r (solid line). A point q £ R(C r ) PI CV is also plotted. 

exhibits a homoclinic transversal intersection if their stable and unstable manifolds intersect 
each other and the tangents of the manifolds are not colinear at the intersection. Transversal 
homoclinic intersections (or homoclinic tangles) beat at the heart of chaos. This is because 
in the region of a homoclinic tangle, initial conditions are subject to a violent stretching and 
folding process, the two essential ingredients for chaotic behavior. 

Let us start to show the main result of the letter considering some important charateristics 
and properties of the flows generated by the velocity fields v and v. The flow <&j associated 
with the autonomous velocity field vo is defined for every (x, y) ^ (0, 0). Thus, given a positive 
time T , the map R° = $^ : 5R 2 \ {(0, 0)} -> 5R 2 \ {(0, 0)} is well defined. It is straightforward 
to show that the map R° written in polar coordinates results R° = (R°(r,9), Rg(r,9j) = 
(r, 9+ -Tq). Map R° keeps invariant the set of points with same radius; i.e., it keeps invariant 
the circles. Moreover, the circles are rotated with a rate of rotation inversely proportional 
to the radius [see Fig. 01a) and (b)] . Observe that map R° resembles a twist map with the 
difference that in the present case, the rate of rotation grows to infinity when the radius is 
reduced. In this respect, it is important to mention that for generic conservative perturbations 
of the twist map the existence of homoclinic points associated to a saddle periodic point was 
proved [20] . This result was also extended to time periodic perturbations of a flow which 
exhibit an elliptic singularity [21]. 

We have assumed that the vortex moves periodically along a curve, that is, r v (t) — r v (t + 
To). Then, the time dependent velocity field v induces a flow $ t which is defined for every 
(x,y) ^ {(0,0)}. So, the map R = $ To : 9? 2 \ {(0,0)} -> 5R 2 \ {(0,0)} is well defined. Note 
that flow $ t is generated by a time periodic Hamiltonian [Eq. (|10fl ], so R is a conservative 
map. Also, R can be extended continuously to (0,0) defining i?(0, 0) = (0, 0) [note that R — > 
when (x, y) — > (0, 0)]. Then, it follows that map R verifies: 

Property A: Like R°, map R also has the property that circles are rotated with a rate of 
rotation inversely proportional to the radius [see Fig.QJa) and (b)]. In other words, if map R 
is written in polar coordinates R(r, 9) = (Rr(r, 9), Rg(r, 9)), then d r Rg is of the order of 1/r 2 . 

Property B: The image of a small circle of radio r intersects transversally this circle; i.e, 
R(C r ) H C r ^ and if q 6 R(C r ) n C r then the tangent to R(C r ) and to C r in q are not 
collinear [see Fig.^c)]. 

From properties A and B of the perturbed map R it follows that for arbitrarily small r 
the map R has a fixed point po with radius smaller than r exhibiting homoclinic transversal 
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Fig. 5 - Schematic plots to show the existence of a saddle periodic point po of the map R with a 
homoclinic transversal intersection, (a) The curve r — * (r,6o(r)) is plotted with thick solid line, and 
the curve r — * R((r,0o(r))) = (r, 6i(r)) is plotted with thick dashed line. Two invariants circles (C ri 
and C T2 ) and their perturbed mapping are also plotted. Note the transversal intersection between 
the circles an their respective mapping (points qi and q2). (b) The rectangle B (shaded area) is 
bounded by two pieces of arcs C\ and C2 contained in C ri and C T2 (where n and ri are closed to 
ro verifying n < ro < n) and two segments ii and h contained in two different rays of constant 
angle. The segments 7„ and y 3 , that connects the periodic point po with the transversal intersection 
qo, are ploted with dash-dotted lines. The mapping R(B) of the considered rectangle is also ploted 
with thick dashed line. 

point. Of course, this result implies that the vector field of Eq. 10 has a saddle periodic orbit 
with an homoclinic transversal intersection. 

We show the previous result in two steps. First, let us display that the perturbed map 
R has a fixed point. Property B guarantees the existence of the curve (r, 9o(r)) plotted in 
Fig. Ufa). Note that 9o(r) is the angular coordinate of a point on the intersection of a circle 
C r with its image R{G r ). A point on such a curve is mapped to other point with the same 
radius r but with a different angle 9i(r), 

R(r,6 (r)) = (Rr(r,6 (r)),R e (r,6 (r)) = (r,0i(r)). 

It is clear that map R has a fixed point if there exist an ro such that ^o(^o) = $i( r o)- I n other 
words, curves (r, 9$(r)) an d ( r ; ^i( r )) °f Fig-EI a ) intersect at ro. From property A follows that 
the variation of #i(r) is larger than the variation of 9o(r); in fact, the derivative of 9\(r) is of 
the order of l/r 2 . This fact guarantees the existence of ro in the vicinity of r — * [19]. 

Now we will see that the mentioned fixed point has a transversal homoclinic intersection 
of their stable and unstable manifolds. We recall that the stable manifold is the set of points 
that converges to the fixed point by forward iteration of the dynamic. Conversely, the points 
on the unstable manifold converge to the fixed point by backward iteration. In Fig. \5j[b) we 
consider a circle Co of radius ro containing the fixed point po- Due to property B, Co and 
its image R(Co) has al least an additional point of intersection denoted by go- Points po and 
qo are connected by the segments 7 S and "f u of Co and R(Co) respectively. Let us consider 
the shaded rectangle B of Fig. OJb) that contains the fixed point po- Using property A, we 
have deduced that segments of constant angle l\ and I2 are stretched by R, and this shows 
that rectangle B is contracted by R along directions close to the tangent of the circles C r and 
expanded along vectors close to the tangent of R(C r ). This is displayed in Fig. OJb) where 
the mapping of the rectangle B is plotted with thick dashed line. Moreover, we have proved 
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in Ref. [19] that R{B) is close to -f u and intersects Co in a point near q , and that, R^ 1 (B) 
is close to 7 S and intersects Co near q . This implies that the unstable manifold of p is close 
to 7„ and the stable manifold is close to 7 S . We recall that -f u cross -f u at q , then the stable 
and unstable manifolds of po have a transversal intersection near qo . 

In summary, we have found an universal mechanism leading to quantum trajectories having 
chaotic behavior. We have shown that the movement of vortices is a generic time dependent 
perturbation of an autonomous velocity field which creates a saddle periodic orbit with a 
transversal homoclinic intersection of their stable and unstable manifolds. This transversal 
intersection generates the well known Smale horseshoe which is the origin of complexity . 
Our results should be useful due to the fact that such deterministic quantum orbits are an 
important theoretical tool for understanding and interpreting several processes in different 
fields. On the other hand, our geometrical analysis of a singular velocity field could be 
important for both theoretical and applied problems of dynamical systems, as for example 
advection in non-stationary fluids [22]. 
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